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Abstract 

Strongly consistent and asymptotic normal estimators of the Hurst index of a stochastic differential equation driven 
by a fractional Brownian motion are proposed. The estimators are based on discrete observations of the underlying 
process. 
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1 Introduction 

Recently, much attention has been paid to the study of almost sure convergence and asymptotic normality of the first 
and second order quadratic variations of a wide class of processes with Gaussian increments (see [?], [?], [?], [?], [?], [?], 
[?], [?]). The conditions for this type of convergence are expressed in terms of covariances of a Gaussian process and 
depend on some parameter 7 € (0, 2). These results help us to obtain the estimate of the parameter 7 and to study its 
asymptotic properties. If the process under consideration is the fractional Brownian motion (fBm) then this parameter 
is the Hurst index. 

The aim of this paper is to obtain a strongly consistent and asymptotically normal Hurst index estimator when the 
solution of a stochastic differential equation (SDE) driven by a fBm with the Hurst parameter 1/2 < H < 1 is observed 
at discrete points. For the first time, to our knowledge, the estimator of the Hurst parameter of a pathwise solution of 
a SDE driven by a fBm and its asymptotic in some special cases were considered in [?]. A more general situation was 
considered in [?]. Estimates for the Hurst parameter were constructed according to first- and second-order quadratic 
variations of the observed values of the solution. Only a strong consistency of these estimates was proved. The rate of 
convergence of these estimates to the true value of a parameter was obtained in [?], when the maximum of lengths of 
the interval partition tends to zero. 

In this paper we offer estimators which are different from the estimators used in previous papers. In general, these 
new estimators allow us to obtain not only strongly consistent but also asymptotically normal estimators. 

The paper is organized in the following way. In section [2] we present the main results of the paper. Section [3] 
is devoted to several known results needed for the proofs. Section [S] contains the proofs. Finally, in section HJ] some 
simulations are given in order to illustrate the obtained result. 



2 Main result 

Consider a stochastic differential equation 

Xt=£+ f f(X s )s+ f g(X.)Bf, te[0,T], T>0, (2.1) 
Jo Jo 
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where / and g are measurable functions, B H = {Bf 1 : ^ t ^ T} is a fractional Brownian motion (fBm) with Hurst 
index 1/2 < H < 1, £ is a random variable. It is well-known that almost all sample paths of -B w have bounded 
p- variations for p > 1/H. 

For < a < 1, C 1+Q (R) denotes the set of all C^functions R -> R such that 

Bup| fl '(s)| + B up ^'^)-9'(y)\ <0o 

Let (fi, J r , P,F), F = (J-t) tg [o,T] be a given filtered probability space, and let £ be J^-measurable r.v. Assume that / 
be a Lipschitz function and g € C 1+a (R), jj — 1 < a < 1. It is well known (see [?], [?], [?], and [?]) that for the fixed 
fBm (B H ,F), 1/2 < H < 1, under the above assumptions, the right-hand side of (|2.ip is well defined and there exists a 
unique F adapted process X — (X t )t£[o-,T] with almost all sample paths in the class of all continuous functions defined 
on [0, T] with bounded p-variation for any p > i and such that 

• X = £ a.s., 

. P f f*|/pr,)|B + I f' g(X,)B?\ <oo)=l for all * G [0; 31. 



We first construct a strongly consistent and asymptotically normal estimate of H if we have additional information 
about function g in (|2.1[1 and observed values of the solution on discrete points. 

Theorem 2.1. Assume that function g is known and there exist a random variable £ such that P(? < oo) = 1 and 

1 

SUP \n(X \\ ^ ? °- S - 

te[o,T] \g(^t)\ 



Then 



- rp , 



H n — > H a.s., 
2^i\n^{H n -H) ^N(0;a 2 ) for H 6 (1/2,1) 

with a known variance a 2 defined in Subsection \3.1\ where 

fl '^[-E[ g(L ; ) j J fom>T, A (2) X± T = X j _ T -2X 1 _ LT + X 1 

/T\ 2x 2 -1 

<p n ,T(x) = y — J (4 — 2 ), (/p n T is t/ie inverse function of(p n> T, x G (0, 1), n > T. 

It is natural to try to find conditions when it is possible to refuse restrictions of the previous theorem. For this purpose 
we need some additional definitions. We assume that the process X is observed at time points -2—T, k = 1, . . . , m n , 
where m n is an increasing sequence of natural numbers defined as follows. Take measurable non-decreasing function 
f : (0; oo) — \ (0; oo) which satisfies condition 

V5 > f{x) = o(x s ), x — > oo. 

For example, one can take take <p(x) = 1 or ip(x) = ln a n, a > 0. Let's assume, it is known that H £ (/?; 1) for some 
P £ (§; 1). Then we take k„ = [n 2 ^(fi(n)]. If we only know that H 6 (1/2, 1) then we set k n = [ni£>(n)]. Put m„ = nk n . 



k 2 ' 

Next denote 



h[t) = 5 2 (X t ), ft = sup h(t); I n = {-k„ + 2j : j = 2, . . . , fc n }, 

tS[0;T] 

fen 



j=-k n + 2 



where 1 ^ A; < n - 1 and tf = ^-T. For each n ^ 1 define A n = { ^T, §T, . . . , ^T}. Finally, set 

1 1 w {1 l 

H ^ k) =2-2W2 1 *^ k = l,...,n-l. (2.3) 

n,fc 

The following theorem holds. 

Theorem 2.2. Assume that P(k > 0) = 1 and a sequence of random variables (r n ) has the following properties: 
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(i) r n G A„ Vn; 

(ii) P(liminfft(T n ) > 0) 
Let in = ^ T n . Then 



2 In 2y/K (H n (i n ) - H) -i+ 7V(0; a%) 



(2.4) 



Hn{in) — ► H a.s. and 

uiit/i a known variance a\ defined in Subsection \3.1\ 
Let us give several remarks on the obtained result. 

Remark 2.1. Note that we impose very mild assumption on the model since condition k > a.s. simply states that 
for almost all sample paths function t h-> g(Xt) is not trivial one. 

Remark 2.2. Considering theorem 12.21 it is natural to pose question how to construct a sequence satisfying (ii). 
Proposition 12.11 gives several possible solutions however one might consider the other alternatives. It is straightforward 
to observe that requirement (ii) is nothing more than a statement that for almost all sample paths of the sequence 
(h(r n )) has to be bounded away from zero. Probably the most simple case is the case when h(t) is a.s. bounded away 
from zero by a non-random constant. Then one can fix any to G (0; T) and define (r n ) by relation 



Wn — to = min 

l^tSn-l [ 



kT 



-to 



Then r n — > to as n — > oo and conditions (i)-(ii) are satisfied. 

Proposition 2.1. For each n > 1 put = E?=i W^f , m = 1,2, and define 



: argmax W^l; 



,•(2) 



i(3) 



- argmm 

l<k<n-l 



argmm 



argmin 

l<fc<n-l 



■ T. 



W, 



(i) 



+ 



(2) 



w. 



(1) 



w 



<2) 



IT 



+ 



'(2) 



(2) 
n,fc 



w 



a) 



IfP(K > 0) = 1 then the sequences of random variables (Tn), j = 1, •••,4, satisfy conditions (i) and (ii). 

Remark 2.3. Note that in usual case the sets argmax and argmin in definition of i„ consist of one integer point 
k G {1,2,..., n}. If this is not the case we assume that i„ is equal to the smallest value of the corresponding arg set. 

3 Preliminaries 



3.1 Several results on fBm 

Let A«£f = Bf - 2^ + Bg,, i = 2, . . . ,n and ^ B H (n) = (A^Bf, A^Bf, A«B|) , i.e. A( 2 >B H (n) 

is a centered Gaussian vector of second order increments. We denote by \i, n ,i = 2,...,n, the eigenvalues of the 
corresponding covariance matrix whereas A* stands for a maximal eigenvalue. Applying standard multivariate theory 
one obtains 



E|| A (2) B H (n)f = ]T Ai,„ = (n - 1) 



4- 2 2 



Var||A (2) B J/ (n)|| 2 = 2^A 2 , n , 



where ||a;|| = \f x\ + ■ ■ ■ + x\, x G R fc , is a standard Euclidian norm. In [?] it was shown (see also [?]) that 



Var||A (2) B H (n 



O 



and A* = O 



(3.1) 



In what follows we make use of quantity 



V n . 



4-2 2H 



Y^(T- H A {2) Bl T )\ 



(3.2) 
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By self-similarity of fBm one deduces that 



„2ff-l 

T7 d 71 II A (2) r,H I sl|2 

K,t = - — ^777 IIA k 'B (n) 



4-2 2ff 

The latter relationship enables us to estimate deviation of V n ,T from 1. To do this we apply inequality given below. 
It follows from (|3.1|l and main result of [?] 



V ( \\\A {2) B H {n)f -E||A (2) B H (n)|| 2 > e) 



^ 2exp 



Cie C 2 £ 2 



for large n. In view of the above we have the bound 



^ 2 exp — min ( Ci 



2H 2 4if-l\l 

ten , 02£ n ) 



n- 1| 



> £ 



|A (2) B H (n)|| 2 -E||A (2) B H (Ti)|| 2 | > 



(4-2 2H )e 



<2exp[-min(Cie(4-2 2H )n,C 2 e 2 (4-2 2i;f ) 2 rc)] sC 2 exp [-C 2 £ 2 (4 - 2 2H )n] , (3.3) 



valid for £ £ [0; J-Al). The other important statements was also proved in [?] (see also [?]). It asserts that 



with a positive definite matrix 



V n , T - 1 

V 2n ,T ~ 1 



JV(0;E H ), 



(3.4) 



(3.5) 



a 2 a 2 /2 

The explicit expression for Eh can be obtained from [?] . It was also given in [?] . To avoid the reader from doing tedious 
calculations and make paper more self-contained we restate the expression given in [?]: 



Pl {l) 



( 2H(2H - 1)(2H - 2)(2H - 3) 
{ 4-2 2ff 



02(H) 



7-3 2 



4_ 2 2» 



12 - 2| 2 " 7 - 41/ - II 2-7 + 6I/I 2 " 7 - 41/ + II 2 " 7 + 1/ + 2 



2-7 



(7-2)( 7 -l)(7 + l) 



■, /ez, 76(0,1); 
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} +c 1 (H)J2p2-2h(1)p2-2h(1-2); a\ = 2y/c 2 {H) + ci(H) ^ P2 -2h(1) P2-2h(/ - 1); 

1=2 1=2 

oo 

a 2 =2 + c 2 {H) + Cl {H) J2 (p2~2h(1)) 2 ; °l = 2~ 2H (3a 2 + aj + 4a 2 2 ) . 



Remark 3.1. Note that application of the delta method to (|3.4|l implies 



Van,a 



JV(0;ffn) 



and 



77| Inl^-lnl 

Vn.T 



A AT(0; crfj) with a 2 H = ^a 2 - 2a 2 . (3.6) 



3.2 Variation 

In this subsection we remind several facts about p-variation. For details we refer the reader to [?]. Recall that 
p-variation on [a; b] of a measurable function / : [a; b] — > R is a quantity 

= sup^|/(x fc ) - /(as fc -i)r, 



where p £ (0; oo) is a fixed number and supremum is taken over all possible partitions 

x : a = xq < x\ < ■ ■ ■ < x n = b, n ^ 1. 

By BV p ([a; b]) denote the set of all measurable functions / : [a; b] — > R having finite p-variation on [a; b]. In the sequel 
we use one well known relationship, namely, 



f,heBV p ([a;b]) 



fh-f(y)[h(b)-h(a)} 



^C p>q V p (f;[a;b])V q (h; [a;b]) 



(3.7) 



inequality is valid for all e > 0; constants Ci, C2 are positive and do not depend on n 
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withC p ,, = C(p- 1 + g- 1 ) = E^ 1 "" ( ^ +q ~ \ | + |>l,andV 7 (s;[o;6]) = (iv(s;[a;6]))r. 

We also make use of the fact that for each e £ (0; 1/2) almost all sample paths of fBm B H = (Bf^ ,t £ [0; T]) belong 
to BVh b ([0; T]), where _ff £ (1/2; 1), H E — ^_ e l and there exists almost surely finite r.v. G £ ,t with the following 
property: 

V H AB H ;[s;t}) ^G e , T \t-s\ H - e , Vs,te[0;T] a.s. 

Note that the last inequality implies the same relationship for the initial process X , i.e. there exists almost surely finite 
r.v. L Ei t for which 

V HE (X;[s;t])^L s , T \t-s\ H -", Vs,t € [0;T] a.s. 



4 Proof of Theorem 12.11 

Before presenting the proof of this theorem, we want to make few comments on the notion of symbols O and o. Let 
(Y n ) be a sequence of r.v. and an increasing sequence of real numbers (h(n)). Assume that there exists a.s. finite r.v. 
C with the property 

\YJ < 



h(n) 



In this case we write Y n = 0(h 1 (n)). If one has an array (Y n ,k), k — 1, . . . , n, n ^ 1, then Y n ^ = 0(h (n)) means 
that there exists a.s. finite r.v. ? satisfying 



max \Y n . k \ ^ —— 
l^fesjn h(n) 



If Y n — <;„ ■ h 1 (n) and c„ — > a.s. asn->oo, we write Y n — o(h 1 (n)). Similarly, if \Y n . k \ =s n -h 1 (n) and 

? n — > a.s. as n — > oo, we write Y n ^ = o(fa (n)). 

Lemma 4.1. Setti — —T. The following (expansion, decomposition, asymptotic) relation 



A^X U =X ti -2X ti _ 1 +X U _ 2 = g(X ti _ 1 )A m B? i +0 



1 



,2(H-e) 



i = 2, . . . n. 



holds 

Proof. Since 



^'X u = f{X s )s- f(Xs)s+ g(X a )B:- 

Jti-l Jt-i-2 Jti-1 

an application of Lipschitz-continuity of / and boundedness of variation of X yields 



i = 2, 



/(*.)? 

■1 

2T# 



f(Xs)s 



\f(X a )-f(X tl _ 1 )\s- 



1/(^.0 -/(X s ) Is 



/2T 

max sup |X S — Xt i _,| < KL e .T I — I , 

n l^is?" t^^ssji, V n / 

where if is a Lipschitz constant and L e ,t is a r.v. bounding variation of X (see Section [3}. 
Next we observe that 



By JT7) 



5 (X s )Bf - / g(X B )B» = / S (X S ) - S (X ti _jBf - 

+ s(x ti _j (ab^-ab^ 



liXJ-giX^JBi 



g(X s ) - giX^Bf 



i = 2, 



;\\g'\\o C Hs ,H s V„ e {X',[U-i',U])Vh. [B H ;\U-rM 



^\\g'\\aoCH € ,H € L Et TG E ,T ( — 



2(H-e) 



(4.1) 



where L Ej t,G £ ,t are a.s. finite r.v. defined in Subsection 13.21 and Ch c ,h c is also finite provided e is sufficiently small. 
The same inequality holds for the term J t * _I [g(X s ) — 3(Xt 4 _ 1 )]B^ r . Hence (|4.ip is true. 
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Proof of Theorem \2.1l Observe first that the function 

' T\ 2x 



p n ,T(x)=(±) X (4-2 2 *), 16(0,1), 



is continuous and strictly decreasing for n > T. Thus it has the inverse function f n 1 T for n > T. By Lemma |4. II we 
obtain 



2 . 



;i) 2H (4-2^)v— 



T) 

2 



E(a (2) 5| t ) 2 + o(^L pt ) 



4- 2 2H 

i=i 



™ 9 1 

£)(r-^flf r ) +0(^). (4.2) 



Note that 



<p n ,T(H n ) , 

v 1 a.s. as n-> oo. 



(p n ,T(H) 

Thus for every <5 > there exist no such that 

<fin,T{H + 5) < <p n , T (H„) < ifin )T (H - S). 

for n ^ no, i.e. 

<p n ,T(H n ) G (tp n ,T(H + S), tp n ,T(H - 5)) Vn > no- 

The function ip n ,T is strictly decreasing for all n > T. Thus H n 6 (i? — <5, + 5) for all n ^ no and it is strongly 
consistence. 

Now we prove asymptotic normality of H n . Note that 



V5n,T(H) 7 VT/ 4-2 



2 if 



2(ff„ - fl) In (J) + In (4 - 2 2H ") - In (4 - 2 2H ). (4.3) 
By Lagrange theorem for the function h(x) = ln(4 — 2 2x ), x £ (0, 1), we obtain 

In (4 - 2 2Hn ) - In (4 - 2 2H ) = (H n - H) [ ti (H + t(H n - H))dt 

Jo 

The equality (|5.1[) we rewrite in a following way 

Vn 1*- ■ 

It follows from (|3.4|l and Remark 13. ll that 



1 = 2 



iV(0;cr 2 ) 

for sufficiently small e. 

5 Proof of Theorem [231 

In what follows we use notions already introduced in Section [2] without restatement. In order to prove our main results 
we need two auxiliary lemmas. Let 

(2k ^j 211 ' 1 k ™ ( h h h \ 2 

V ; j=-fc„+2 



2H-1 



with 1 < k < n - 1. 
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Lemma 5.1. The following relations hold 



max 

l^fe^n-l 



a.s. 



v$(fc) 



N(0;a H ), for each k = 1, . . . , n — 1; 



max \V^ p i(k) - 1| — ► a.s. /orp = 1,2; 

l<fc<n-l 



max 

S;,le{l,...n-1} 



a.s. /or p = 1,2; 



where a% is the same as in (13.611 . 

Proof. We will give the proofs of the first two assertions. The proof of the rest is obtained by obvious alterations. 
By self-similarity and stationarity of first order increments 



v$(fc) 



„2H (off -9B h 



B H 



3=2 I ^3 



V I B H — 2B H + B H I o- 2H ^P kn I B H ~2B H 4- B H 



2 2H ~ 1 T4„,i 



with V"2fe„,i, Vfe„,i defined by (|3.2fl . Hence 



Let Cn t = max 

l<fc<n-l 



Then for all e > 



P(Cn,T > e) < nP 



V 2h „ 



-1 >e 



Let (e n ) be a sequence of positive numbers such that e„ J, 0, n — >• oo, and let 2fe n 1 < e n . Since — 1| ^ 2e n for all 
x,y £ (1 — 2e„; 1 + 2e n ) application of (|3.3|) gives 



- 1 



■ • 2.-,, | 1' ( > -;, ... ~ J| 2r„ ) +I>( |Vfc„,i - J| ■ 2.:-,, 

P 



t r k n 1 

Vfc_,i ; > £r 



<4exp (-Ce n fc„) 



Set g n = {s/Ckn ) x \/3 In n. Then 2fc„ 1 < e n for large n and the series X^P(C™,t > £n) converges. Consequently, 
Cn,T — s* a.s. 

We have already proved that the distribution functions of 



1 and ykn 



V2k n ,l 
Vk n A 



are equal. By Remark 13. II we obtain the required assertion. 
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Lemma 5.2. If e € (0;0.3) is sufficiently small then the following relations hold: 



(1) _ / k \ T"(4-2") (1) / 1 \ _ / k T \ g««(4-2" 

^ In j (2fe„)"-i ^• TW + U H - e ^ H -V" U . 



(2) _ / fc \ T"(4-2") [2) ( 1 \ _ / fc \ r"(4-2"*) 



.2H-1 ' 



uiit/i 1 ^ fc ^ n — 1. 

Proof. Step 1. In the same way as in Lemma [4. II one obtains relation 

A (2) X tr = Xty. - 2X tT _ l + X tT _ 2 = g(X tT i )A (2) BtL + O ( ^(g-.) ) > * = 2, . . .m». 
Step Note that 

|ft(t) - =|fl(X t ) - g(X s )|| S (X t ) + 5 (X S )| 

^2[| S '| 00 VH e (X; [0,T]) + | fl (0|] IffVl^t - *.| 
apff'IccV^CX; [0,T]) + | fl (0|] |ff'|oo^«.r|* - «| K ~", 
where |(/|oo = su P tg [ .T] • ^ follows from relation (|5.5p . inequality (|5.6[1 . and Lemma l5.1l that 

<i= e (A (2 »x tr+ , T ) 2 = e ^fc+^jfA^jVo^ 

j=-fc„+2 j=-fc„+2 ^ ' 



k " 



+ o 



J=-fc«+2 j=-fc„+2 



2H/,i o2H\ / /ot \2H-1 fc « 



n J (2fc„)^-i I (4-2") ' 

+ o( fc " ^ + o ^ fcn 



H 

, + k T 



=h[ — T ' 



/* \ r 2g (4-2^) (1) / i \ 



Reasoning as in the previous case, we can prove that 

n2HfA o2-H^ 



Next we consider the quantity W„ . Note that 



V — — !— h ( — t\ r"(4-2") (1) / 1_ 



rp2H _ 22H 



(2fc„) 2H 

+ o 



Since almost all sample paths of X are continuous, it follows that 

/ ;„ \ fT 



Since 



T 



^ max V, 



(i) 



(5.8) 



then by Lemma 15. II and (|5.T|) we get that the right side of (|5.8|) converges to 0. Thus 



W (i) = T 2H - 1 (4-2 2 



T 



and relation (|5.3[1 holds. (|5.4|l follows similarly. The proof of the lemma is thus complete. 

Proof of Theorem HO Fix u G A = {n > 0} n {liminf ft(r n ) > 0}. Next find 5 > and n < oo for which ft(r n ) > 5 
whenever n no. It follows from Lemma |5. II that 

r(p) 



V^(in) >1, p=l,2. 



This, together with (|5.1|l and (|5.2|l . yields 



InW, 



ci) 



lnV; (1 ^(i„) - (2i7- 1) In 2 + In 



rp2H _ 22-ff \" 



h (t„) 



and 



InW^ = InV$,(i„)+ln 



(2), 



^ (T n ) 



rj-i2.II q2.H 1 



+ In 1 + O 



+ In 1 + O 



(fen) 211 " 1 

By applying Maclaurin expansion, we obtain that In (l + 0( n u- e )) = Q( — ?r=i )■ Thus 



H " ( * n) = I ~ 2^2 ( lnW/ ^^) " = H ~ ^ 

Using Lemma |5. II we have 



■In 



ln2 VV«(in) 



+ o 



(5.9) 



Since P(j4) = 1, the consistency of H n (i n ) is proved. 
Now rewrite equality ()5.9|) in the following way 

n-, , s /t- /^tWA / 

21n2v^(ff«(i™)--ff) = -VfcTln +° V « 

For sufficiently small e the last term goes to a.s. This, together with Lemma |5. II and the delta method (see Remark 
13. ip . yields the required statement. 

Lemma 5.3. Set A — {k > 0} n {X ts continuous} and let 



->• 1 a.s. 



ro(w) 



To(w) = 



0, for U! G A c , 
inf {t G (0, T) : ft(t 

0, /or w G A c 

inf {t G (0, T) : Th(t, oj) = J Q T /i(s, w) ds}, /or weA 



(5.10) 



T/ien to and to are T -measurable r. v's. 

Proof First we prove that ro is measurable r. v. The proof follows from the equality 

{u : t (u) ^ x} = ( Pi (J |cj G ^4 : j/i(g,oj) — sup /i(s, < — | j U A c . (5.11) 

Now we verify it. Denote the left side of (|5.11[) by B and the right side by C. First we prove that B C C. If uj G B n A c 
then it is obvious that w 6 C. If w 6 B (1 A then h(jo{uj), uj) = sup s<T h(s,uj) and there exists a sequence of rational 
numbers (q n ) G Q fl [0, x] such that q n — > To(ui) and 

|/i(g n ,w) — Ii(to(uj),lu)\ < — for all n > 1. 
1 1 n 



Thus uj £ C and B C C. 
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Now we prove C C B. Assume that lo £ C n A. Then for every n> 1 we can choose g„ £ Q Pi [0, x] such that 

\h(q n ,uj) — suph(s,w)\ < i . 



The sequence (q n ) is bounded. Thus we can choose a subsequence (q nk ) which converges to some limit xq. It is clear 
that xq G [0, x]. By definition of to(oj) we get that ro(w) ^ xo- Thus u> G B and C C B. 
.F-measurability of r. v. ro follows from equality 



{u) : r (cj) ^ x} = f P| (J jw G A: |%,6j) - f h(s,uj)ds ^ -\^\ U A c 



- ijeQn[0,a: 

and Lagrange's mean value theorem. 

Proof of vrovosition [2A\ For each n determine r n G A n = { — T, 



-T} from equality 



(5.12) 



I kT 

|r«(w) — To(u)\ = min T o{u)\ 

l^fc^n — 1 I n 



(5.13) 



Then t„ > ro and t„ is ^-measurable r. v. for every n. By definition of z„ and (|5.1 



(,<•>) 



rp2H ^ _ : 



+ o 



T 2ff (4 - 2 2 
(2fc n ) 2H ~ 



+ o 



Multiplying both sides by T ^"]_ 22H j and rearranging terms gives h (j'n^J h (r n ) + o(l). Consequently, 

liminf h (t^A > liminf {h (r n ) + o(l)) = re> 0. 



With the remaining sequences (rn ) we proceed in a similar way. Let ro be a random variable as defined in (|5.10|l . 
Next we denote \\h\\ = J Q T /i(t)t. Condition k > implies \\h\\ > 0. A sequence of random variables (r n ) is defined by 
using relation (|5.13|) and replacing ro by rb. By using Lemma 15.21 positivity of \\h\\, and finiteness of sup i<T h(t), we 
obtain 



Wi p) k _Tfc(&T)+o(l) _ Th (~T) 



wi p) 



W 



(2) 
n 

(2) 



w 
w 



(1) 



IF 



\\h\\+o(l) \\h\\ 

Kl _(il 2H ^2MM±°( 1 ) 



| +0(1) 



o(l), p = l,2; 
"i-i Th (k T \ 



Th($T)+o{l) „ H _,Th(±T) 



(1) /1\2H--1 



(i)^" 1 ||/l||+ 0(1) 



HI 



+ o(l); 



IF 



■(2) 



0(1). 



Therefore by definition of i„ and r. 





- 1 










^ (1) (3, 


- 1 


+ 









- 1 



Th (t„ 



+ o(l) 



IF 



< 2 



(2) 



+ 



IF 



71, l« 



(i) 



IFA 



(2) 



(2) 



1 =2 



IF 



(i) 



IF 



Th {T n )-\\h\\ 



I h\\ 



+ o(l) 



(2) 



+ 1 2 



T/i (r„)-||/i|| 



/ill 



IF 



+o(i: 



+ 



if 



<1) 



Since h is a continuous process then convergence r n — > ro implies Th(r n ) — > Th(ro) = \\h\\. Consequently, right-hand 
side of each inequality vanishes. So does the left-hand side and it remains to observe that in such case lim n ->oa h{jn ) = 
j = 2, 3, 4. 
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6 Simulations 



The goal of this section is to compare the performance of the obtained estimators. The processes chosen for simulations 
were 

I: X t = X + f* sin(X t ) dt + J* cos(X t ) dB? , t G [0, 1], X = 1; 

II: Xt = X + f* sm{X t ) dt + f* (2 + cos(X t )) dflf , t G [0, 1], X = 1. 

The sample paths of the fractional Brownian motion were generated using the circulant matrix embedding (or Wood- 
Chan) method (see, f. e., [?]). The sample paths of the processes defined above were obtained by applying the Milstein 
time-discrete approximation. The sample path lengths and the values of the Hurst index were chosen to be, respectively, 
n 2 + l where n G {50, 100, 150, 250, 500} and H G {0.55, 0.65, 0.85, 0.95}. For each tuple (n 2 + l, H) we simulated 200 
sample paths of both considered processes. The estimators H™ , i G {1,2,3,4} were calculated using the formula (|2.3|) 
with the sequences i n chosen as described in the Proposition 2.1 and were calculated for the process (I). The estimator 
H$ corresponds the statistic of the Theorem 12.11 and, due to the requirement for the function g(Xt) to be separated 
from zero, it was initially calculated for the process (II). However, after comparing the numerical characteristics of H£ 
to those of H", i G {1,2,3,4} we decided to check how it performs in case the requirement for the function g(X t ) is 
not fulfilled, hence we calculated H$ for the process (I) as well. For notational simplicity, let Hi := H™. 

The table [1] presents the mean squared errors and the average absolute deviations of the considered estimators 
(* - MSE was multiplied by 10 5 instead of 10 3 ). It can be seen that the numeric characteristics of H3 and HA are 
considerably better than those of HI and H2. The estimator HI displayed a rather large and slowly decreasing negative 
bias. Additionally, H3 and HA show similar behavior for H < 3/4; in fact, on the average for 80% of the sample paths 
the estimated values of H3 and HA were identical. However for larger values of H the estimator if 3 was noticeably more 
stable. The numerical characteristics of the estimator H5 surpassed those of HI — HA for both considered processes, 
yet in case of (I) there were multiple outliers. The magnitudes of these outliers were directly proportional to the 
lengths of the sample paths. These observations are further illustrated by figures ?? and ?? presenting, respectively, 
the comparison of the estimators HI, H3 and H5 (for both (I) and (II)) for the fixed values of H and for the fixed 
sample path lengths n 2 + 1. 







MSE(Hi) ■ 


10 3 






E\Hi-H\ 


• 10 






n 


H 


0.55 


0.65 


0.85 


0.95 


n 


H 


0.55 


0.65 


0.85 


0.95 




HI 


33.21 


36.35 


26.61 


22.93 




HI 


1.439 


1.529 


1.262 


1.151 




H2 


33.38 


30.77 


26.10 


19.64 




H2 


1.428 


1.403 


1.238 


1.133 


50 


H3 


9.127 


9.303 


7.228 


6.244 


50 


H3 


0.748 


0.778 


0.670 


0.645 




H4 


8.983 


10.59 


8.874 


9.117 




H4 


0.750 


0.809 


0.744 


0.771 




H5(I) 


1.863 


0.900 


0.379 


4.831 




H5(I) 


0.074 


0.068 


0.054 


0.094 


* 


H5(II) 


0.419 


0.390 


0.207 


0.067 




H5(II) 


0.017 


0.016 


0.012 


0.007 




HI 


13.60 


11.34 


6.307 


7.395 




HI 


0.926 


0.854 


0.625 


0.696 




H2 


12.38 


9.380 


6.057 


4.894 




H2 


0.874 


0.810 


0.624 


0.554 


150 


H3 


2.506 


2.047 


1.418 


1.043 


150 


H3 


0.395 


0.340 


0.299 


0.259 




H4 


2.594 


2.216 


2.013 


1.717 




H4 


0.401 


0.351 


0.355 


0.322 




H5(I) 


0.009 


0.024 


0.159 


0.042 




H5(I) 


0.007 


0.010 


0.019 


0.011 


* 


H5(II) 


0.030 


0.027 


0.017 


0.006 




H5(II) 


0.004 


0.004 


0.003 


0.002 




HI 


4.128 


3.337 


2.941 


1.970 




HI 


0.525 


0.463 


0.435 


0.357 




H2 


3.176 


2.581 


1.721 


1.573 




H2 


0.454 


0.400 


0.334 


0.318 


500 


H3 


0.372 


0.334 


0.232 


0.180 


500 


H3 


0.148 


0.138 


0.118 


0.108 




H4 


0.382 


0.365 


0.278 


0.290 




H4 


0.151 


0.145 


0.129 


0.132 




H5(I) 


0.000 


0.003 


0.015 


0.029 




H5(I) 


0.003 


0.004 


0.008 


0.010 


* 


H5(II) 


0.008 


0.007 


0.005 


0.002 




H5(II) 


0.002 


0.002 


0.002 


0.001 



Table 1: Comparison of the estimators 



11 




~i 1 1 r 

0.55 0.65 0.85 0.95 




~i 1 1 r 

0.55 0.65 0.85 0.95 





1 — i — ' ' o 




1 ~*~ 

o 

o ° 



~i 1 1 r 

0.55 0.65 0.85 0.95 



(a) HI, n=50 



(b) HI, n=150 



(c) HI, n=500 




"i 1 1 r 

0.55 0.65 0.85 0.95 



-9- 



'• | -■- ~o~ 

t- 



"1 1 1 1 — 

0.55 0.65 0.85 0.95 




"i 1 1 r 

0.55 0.65 0.85 0.95 



(d) H3, n=50 



(e) H3, n=150 



(f) H3, n=500 



bBB 



— I 1 1 1 — 

0.55 0.65 0.85 0.95 




~i 1 1 r 

0.55 0.65 0.85 0.95 




~i 1 1 r 

0.55 0.65 0.85 0.95 



(g) H5(II), n=50 



(h) H5(II), n=150 



(i) H5(II), n=500 



— I 1 — 

0.55 0.65 



— I 1 

0.85 0.95 



— I 1 1 1 — 

0.55 0.65 0.85 0.95 



1 1 1 1 

0.55 0.65 0.85 0.95 



(j) H5(I), n=50 



(k) H5(I), n=150 



(1) H5(I), n=500 



Figure 1: the dependance of (Hi — H) on the value of H 
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50 100 150 250 500 



50 100 150 250 500 



50 100 150 250 500 



50 100 150 250 500 



(a) H3, H=0.55 



(b) H3, H=0.65 



(c) H3, H=0.85 



(d) H3, H=0.95 



— H-B-i 



~~ i — i — i — i — i — 

50 100 150 250 500 

(e) H4, H=0.55 



-tr -& 



— I 1 1 1 1 — 

50 100 150 250 500 

(i) H5(II), H=0.55 




~i 1 1 1 r 

50 100 150 250 500 

(m) H5(I), H=0.55 




"i 1 1 1 r 

50 100 150 250 500 
(f) H4, H=0.65 




"i r 



"i r 



50 100 150 250 500 
(j) H5(II), H=0.65 




~i 1 1 1 r 

50 100 150 250 500 



-ir 



~~ I 1 1 1 1 — 

50 100 150 250 500 

(g) H4, H=0.85 




"i 1 1 1 r 

50 100 150 250 500 

(k) H5(II), H=0.85 




I 1 1 1 1 — 

50 100 150 250 500 



-Hf' 



-& 



t 



t 



~i 1 1 1 1 — 

50 100 150 250 500 

(h) H4, H=0.95 




"i 1 1 1 r 

50 100 150 250 500 

(1) H5(II), H=0.95 




I 1 1 1 1 — 

50 100 150 250 500 



(n) H5(I), H=0.65 (o) H5(I), H=0.85 (p) H5(I), H=0.95 

Figure 2: the dependance of (Hi — H) on the length of the sample paths n 2 + 1 
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